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Verification of g-factors for lead monofluoride ground state, PbF 

L.V. Skripnikov, 1, 2, 0 A.N. Petrov, 1,2 ’0 A.V. Titov, 1, 2, @ R.J. 

Mawhorter, 3 A.L. Baum, 3 T.J. Sears, 4 and J.-U. Grabow 5 

1 National Research Centre “Kurchatov Institute” B.P. Konstantinov Petersburg 
Nuclear Physics Institute, Gatchina, Leningrad district 188300, Russia 
2 Dept, of Theoretical Physics, St.Petersburg State University, 198504, Russia 
3 Department of Physics and Astronomy, Pomona College, Claremont, California 91711-6359, USA 
4 Chemistry Department, Brookhaven National Laboratory, Upton, New York 11973-5000, USA 
5 Gottfried- Wilhelm-Leibniz-Universitt, Institut fr Physikalische Chemie 
and Elektrochemie, Lehrgebiet A, D-30167 Hannover, Germany 

We report the results of our theoretical study and analysis of earlier experimental data for the 
g-factor tensor components of the ground 2 n lj / 2 state of free PbF radical. The values obtained both 
within the relativistic coupled-cluster method combined with the generalized relativistic effective 
core potential approach and with our fit of the experimental data from [R.J. Mawhorter, B.S. 
Murphy, A.L. Baum, T.J. Sears, T. Yang, P.M. Rupasinghe, C.P. McRaven, N.E. Shafer-Ray, L.D. 
Alphei, J.-U. Grabow, Phys. Rev. A 84, 022508 (2011); A. Baum, B.S. thesis, Pomona College, 
2011]. The obtained results agree very well with each other but contradict the previous fit performed 
in the cited works. Our final prediction for g-factors is Gm = 0.081(5), G± = —0.27(1). 


INTRODUCTION 


MOLECULAR HAMILTONIAN 


Lead monofluodide, PbF, molecule is one of prospec¬ 
tive systems to search for the electron electric dipole mo¬ 
ment (eEDM). It was studied and discussed during three 
decades in many papers including (THi]- It was recently 
shown in Ref. [7j that some “enhanced” (coincidental) 
near-degeneracy for the levels of opposite parity in the 
ground rotational state J = 1/2 for 207 PbF of the ground 
electronic state 2 n 1 / 2 3 takes place that is caused by 
the near cancellation between the shifts in the energies 
of these levels due to omega-type doubling and the mag¬ 
netic hyperfine interaction. This can lead to suppression 
of systematic errors in an experiment. 

In Ref. Q we have calculated the parameters (more 
generally, the characteristics of atoms in compounds @- 
|n|]) required to interpret the experimental energy shift 
in terms of the eEDM and other effects of simultaneous 
violation of space parity (P) and/or time-reversal invari¬ 
ance (T) including the P-odd anapole moment Q and 
the T,P-odd pseudoscalar-scalar electron-nucleus neutral 
current interaction for the ground 2 Hi/ 2 state. For in¬ 
stance, the effective electric field in PbF was found to 
be greater than or equal to those in the other transition 
element compounds considered (1.7 times larger than in 
HfF + 12, p|, 1.4 larger than in PtH + 141, and 1.1 larger 
than in WC |l5[ and TaN 0))- 

In the present paper our aim is to study the PbF g- 
factor for the 2 n 1 / 2 term which is required for prepara¬ 
tion of experiments on the molecule H, 0, 0 ■ Up to 
now the g-factors have been measured in Ref. 0, 0 
only. Previous theoretical estimations and calculations 
of g-factors have been performed in Refs. QiH- 


We represent the molecular Hamiltonian for 208 PbF 
as 

H mo i = H rot + Hhf s + Hu + H ext . (1) 

Here H rot is the rotational Hamiltonian and Hhf s is the 
hyperfine interaction between electrons and nuclei. Hi 
includes the nuclear spin rotational interaction and also 
effectively takes into account the rotational and hyperfine 
interactions between 2 Hi/ 2 and other electronic states. 
H ext describes the interaction of the molecule with an 
external magnetic field B. Parameters for H rot , Hhf s , 
and Hi are taken from Ref. [uj. For H, : „. t we have: 

H ext = hb B • G ■ S' — (j\ hn B • Ii (2) 

Here S' is effective spin defined by the following equa¬ 
tions: S/|fi >= D|D >, S' ± |D = =pl/2 >= |fi = ±1/2 >, 
Sj_|f2 = ±l/2>=0 (T, H, Ii is the angular-momentum 
operator of the fluorine nuclei, p-B and g-N are Bohr and 
nuclear magnetons respectively, and g i = 5.25773 is the 
19 F nuclear g— factor. 

In the molecular frame coordinate system the tensor 
contractions 

B • G ■ S' = G||B 0 Sq - Gj_(BiS , _ 1 + B_ x Si) (3) 
are determined by the body-fixed g —factors 

G„ =2<* 2 n 1/2 |^- 5 s^|^n 1/2 ), (4) 

Gj_ = {* mi/a \L B + -gs§l\*in_ 1/a ), (5) 

where L e and S e are the electronic orbital and electronic 
spin momentum operators, respectively; gs = —2.0023 is 
a free—electron g-factor; h is the unit vector along the 
molecular axis directed from Pb to F. 
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In this paper the parameters Gy and G± are obtained 
(i) by using Eqs. m from calculation of the electronic 
wavefunction II/ 2 n and (ii) by fitting the experimen¬ 
tally observed transitions reported in Ref. [20]. 


METHODS 


The matrix elements CTol) were calculated using the 
computational scheme similar to that used by us in 
Ref. 0- The basis set for Pb was taken from Ref. Q. 
For F the aug-ccpVQZ basis set [2l[ with two removed 
g-type basis functions was employed. The Pb—F inter- 
nuclear distance was set to 3.9 a.u., which is close to the 
experimental datum 3.8881(4) a.u. [2sj], which was later 
confirmed by Ref. [HI]. Inner core Is—4/ electrons of lead 
were excluded from the correlation calculation using the 
“valence” semi-local version of the generalized relativistic 
effective core potential (GRECP) approach [H|,[25[. Note 
that the approach allows one to account for the Breit in¬ 
teraction very effectively [3 SB, [§7j]. All the other 31 
electrons were included into the calculation. Electron 
correlation effects were considered within the relativistic 
two-component coupled-cluster approach with account¬ 
ing for single and double cluster amplitudes, CCSD, as 
well as single, double and perturbative triple cluster am¬ 
plitudes, CCSD(T). Note that the matrix element ([5]) 
is off-diagonal. Therefore, it was calculated within the 
linear-response two-component coupled-cluster method 
with single and double cluster amplitudes [3]. The 
coupled-cluster calculations were performed using the 
DIRAC12 {29] and MRCC j30] codes. Matrix elements of 
the operators corresponding to m over the molecular 
spinors were calculated with the code developed in Refs. 
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31 34|. 


To obtain the experimental values for Gy and G 1 w e 
have performed two fits using the data from Ref. [20l |. 
In “fit 1” the Zeeman shifts of J = 1/2 to J = 3/2 
transitions for the ground vibrational level of 2 n 1 / 2 elec¬ 
tronic state are obtained by numerical diagonalization 
of the molecular Hamiltonian (H mo i) on the basis set 
of the electronic-rotational wavefunctions. The scheme 
of the calculation is similar to that employed in Refs. 
BE Hi. Only the Gy and Gj_ parameters were opti¬ 


mized. The other parameters of H mo i were taken from 
Ref. [s]. In “fit 2” we have reproduced the scheme de¬ 
scribed in Ref. El. 


RESULTS AND DISCUSSION 

The results of our calculations of g-factors for the PbF 
ground state together with the results of previous stud¬ 
ies are given in Table HI One can see that the value of 
Gy is stable with respect to improvement of the electron 
correlation treatment in the present study (from CCSD 


to CCSD(T) level). 


TABLE I. Calculated values of g-factors (Gy, Gi) of the 2 n 
state of PbF. 


Method 

Gy 

G x 

SCF yj 

0.034 < Gy < 0.114 

-0.438 < Gx < -0.269 

SCF “, [2] 

0.114 

-0.438 

13e-SODCI b , [5] 

0.082 

-0.319 

31e-CCSD, this work 

0.081 

-0.274 

31e-CCSD(T), 
this work 

0.081 


Experiment, [191 

0.12 

-0.38 

Experiment + fit 1, 
this work 

0.081 

-0.269 

Experiment + fit 2, 
this work 

0.085 

-0.271 


“ SCF, self consistent field. 
b 13-electron SODCI, spin-orbit direct configuration 
interaction, :5|. Outer-core electrons 5s 2 5p 6 5d 10 of Pb are 
excluded from the correlation treatment. 


Gy and G± values obtained by fit 1 and fit 2 (see Meth¬ 
ods section) are also given in Table [T] The deviations of 
our fits from the observed Zeeman shifts are given in Ta¬ 
ble HD For the last seven transitions the shifts are repro¬ 
duced with deviations which are much larger than the de¬ 
clared experimental accuracy. One is inclined to suspect 
that the accuracy is overestimated for these transitions. 
We note however, that the experimental (AU/B) 0 \, s val¬ 
ues for Zeeman components that only differ (model in¬ 
dependent) in sign (e.g. Fx, MFx A Fy, MFy = 1, 0 
-A 2, 1 vs. 1, 0 A 2, -1; 1, 1 A 2, 2 vs. 1, -1 A 2, -2; 
1, 1 A 2, 0 vs. 1, -1 A 2, 0; 0, 0 A 1, 1 vs. 0, 0 A 1, 
-1) agree within their error bars, which indicates correct 
accuracy estimations. It is also the case that the devi¬ 
ations for those pairs are systematic and not statistical. 
It seems that the Fx aF(/ = 1a2 pattern is predicted 
to be somewhat too narrow while the Fj, aF;j = 0a 
1 pattern is somewhat too wide. 

We also note that the Gy = 0.085, G_l = -0.271 param¬ 
eters obtained in fit 2 differ substantially from the Gy = 
0.12, G_l = -0.38 values obtained by the same method 
and reported in the Ref. 0. Our results here show 
good agreement between Gy and Gj_ obtained in fit 1, 
fit 2, and the ab initio calculation. While both Ref. flfl ] 
values are higher by a common factor of ~1.4-1.5, the 
origin of the discrepancies is not clear at present and will 
require further investigation. 

Our final values for the g-factors are Gy = 0.081(5) and 
G_l = -0.27(1). It should be noted that these smaller g- 
factor values and their improved accuracy together favor 
the experimental search for the electron electric dipole 
moment and other parity-violating and related effects 
SUzl m PbF due to the additional suppression of sys¬ 


tematic errors. 
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TABLE II. Observed Zeeman shifts (AU/B) 0 b s (MHz/Gauss) 
of the J = 1/2 to J = 3/2 transitions for 208 Pb 19 F [2(|. The 
number in parenthesis gives two standard deviation error of 
the final digits of precision. The subscripts U and L refer 
to the upper and lower energy level of the transition, respec¬ 
tively. F is the total angular momentum of PbF, MF is its 
projection on laboratory axis. The deviation of n-th fit is 
given by 5„ = ( AU/B)bt — (AU/B) 0 ba in units of the last 
digit of precision. 


Unsplit line (MHz) F l 

Fu 

MFz, 

MFc/ (AU/B) ohs [20] 

<5i 

18414.588 

1 

2 

-1 

-1 

0.0665(13) 

-40 




0 

0 

-0.00050(93) 

107 




1 

1 

-0.0635(13) 

9 

18462.193 

0 

1 

0 

0 

0.00032(90) 

-89 




0 

-1 

-0.1369(30) 

-4 




0 

1 

0.1363(29) 

-2 

18497.136 

1 

1 

-1 

-1 

0.00766(86) 

-150 




1 

1 

-0.00729(92) 

119 




1 

0 

-0.1428(17) 

0 




0 

-1 

-0.1328(21) 

-33 




0 

1 

0.1345(13) 

29 




-1 

0 

0.1427(9) 

2 

22574.934 

1 

2 

-1 

-1 

-0.03864(27) 

28 




0 

0 

-0.00005(90) 

2 




1 

1 

0.03851(9) 

-14 




1 

0 

0.1023(36) 

64 




0 

-1 

0.07296(25) 

-264 




-1 

-2 

0.03411(60) 

-211 




1 

2 

-0.03406(49) 

206 




0 

1 

-0.07267(52) 

229 




-1 

0 

-0.10323(79) 

-546 

22691.931 

0 

1 

0 

-1 

0.11114(48) 

371 




0 

1 

-0.11133(41) 

-346 
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